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ABSTRACT 


This  dissertation  presents  a  method  for  reducing  the  number  of  paraunetric  un¬ 
certainties  used  in  the  design  of  a  robust  Hoc  controller.  The  resulting  controller  is 
shown  to  meet  robust  stability  and  performance  requirements  in  the  presence  of  all 
the  modeled  uncertainties.  The  approach  used  involves  grouping  parametric  vari¬ 
ations  affecting  the  same  open  loop  eigenvalue,  then  scaling  one  or  more  of  these 
variations  to  accommodate  the  eigenvalue  chamge  caused  by  all  the  parametric  un¬ 
certainties.  This  method  is  effective  in  those  cases  where  a  large  number  of  parametric 
uncertainties  cause  current  computer  aided  design  software  to  fail  to  find  a  robust 
controller  for  the  plant. 
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I.  INTRODUCTION 


A.  PROBLEM  STATEMENT 

The  advent  of  Hoo  synthesis  and  associated  robust  stability  analysis  techniques 
has  provided  a  process  for  designing  controllers  which  are  capable  of  meeting  es¬ 
tablished  performance  criteria  and  remaining  stable  over  a  wide  range  of  operating 
conditions  and  parameter  uncertainties.  To  be  successful,  however,  this  process  re¬ 
quires  that  all  uncertainties  and  variations  in  the  plant  be  quantified  and  included  in 
the  plant  model.  Modeling  the  plant  in  this  manner  may  result  in  either  a  high  order 
system  or  a  system  severely  constrained  by  the  number  of  uncertainties.  Current 
computer  aided  design  (CAD)  software  often  cannot  find  a  solution  when  the  system 
order  becomes  quite  large  or  when  a  large  number  of  uncertainties  are  present  in  the 
plant  model.  This  does  not  imply  that  an  appropriate  controller  does  not  exist,  only 
that  current  techniques  are  inadequate  in  these  cases.  A  new  technique  is  required 
that  provides  for  the  use  of  current  synthesis  methods  emd  software,  but  results 
in  a  controller  design  which  is  stable  and  hais  acceptable  performance  over  all  plant 
variations.  In  this  dissertation  we  will  present  a  technique  for  finding  a  robust  Hoo 
controller  in  specific  cases  where  the  present  design  methods  fail. 

B.  HISTORY 

The  period  of  modern  control  is  delineated  by  state-space  analysis  of  systems  and 
the  development  of  optimal  Linear  Quadratic  Gaussian  (LQG)  controllers.  The  LQG 
design  techniques,  however,  largely  ignored  uncertainties  in  the  system  and  the  stabil¬ 
ity  margins  of  Linear  Quadratic  Regulator  (LQR)  controllers  with  full  state  feedback 
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proved  elusive  when  combined  with  K&Iman  Filter  estimators  (Doyle,  1978).  Re¬ 
cently,  two  new  techniques  for  controller  design  and  stability  analysis  of  multi- input 
multi-output  systems  have  been  introduced.  Structured  Singular  Value  analysis  and 
Hoo  synthesis. 

Structured  Singular  Value  (SSV)  analysis,  also  referred  to  as  p  analysis,  pro¬ 
vides  a  mechanism  for  determining  the  stability  and  performance  of  a  system  in  the 
presence  of  the  uncertainties  defined  in  the  plant  model.  The  SSV  and  an  equiva¬ 
lent  measure,  the  Multivariable  Stability  Margin,  were  introduced  in  1982  (Doyle, 
1982)(Safonov,  1982).  While  providing  a  definition  for  determining  the  stability  of 
a  system  with  uncertainties,  no  closed-form  solution  exists  for  finding  the  SSV  or 
Multivariable  Stability  Margin.  It  is  possible,  however,  to  calculate  reaisonably  tight 
bounds  on  the  measure  of  stability  margins.  A  considerable  amount  of  effort  has 
been  placed  on  reducing  the  conservative  nature  of  the  upper  bound,  particularly 
since  the  method  for  calculating  the  bound  does  not  distinguish  between  real  and 
complex  uncertainties.  Methods  for  calculating  the  actual  SSV  for  real  variation 
in  the  plant  parameters  using  interative  computational  techniques  (de  Gaston  and 
Safonov,  1988)  and  polynomial  methods  (Elgersma,  Freudenberg  and  Morton,  1992) 
exist,  but  are  unwieldly  auid  not  suited  for  most  design  work.  Recent  research  has 
also  suggested  a  method  for  obtadning  an  uncertainty  model  which  ensures  that  the 
number  of  parametric  uncertainties  included  in  the  model  is  a  minimal  set  (Belcastro 
and  Chang,  1992). 

Hoo  synthesis  provides  a  method  for  the  design  of  a  controller  which  minimizes 
the  peak  magnitude  of  the  closed-loop  frequency  response  of  appropriately  modeled 
systems.  The  first  Hoo  controller  design  methods  were  computationally  difficult  and 
resulted  in  controllers  having  many  more  states  than  the  plant  (Zames,  1981), (Fran¬ 
cis,  1987).  A  breakthrough  in  1988  provided  a  design  technique  which,  similar  to 
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LQG,  involves  solving  two  Riccati  equations  and  results  in  a  controller  of  the  same 
order  as  the  plant  (Doyle,  Glover,  Khargonecker  and  Fr^lncis,  1988).  By  appropriately 
defining  uncertainties  in  the  model  and  adding  performance  weighting  functions  to 
the  system,  it  is  possible  to  design  a  controller  using  Hoo  synthesis  which,  when  tested 
using  SSV  anaJysis  can  be  shown  to  be  stable  and  meet  perform2mce  criteria  over 
all  variations  quantified  in  the  model.  The  controller  design  may  fail  as  mentioned 
above,  however,  when  the  plant  model  is  of  high  order  or  there  are  many  uncer¬ 
tainties  in  the  plant.  Despite  the  recent  advances  in  uncertainty  analysis  and  Hoo 
synthesis,  the  problem  of  finding  a  controller  which  is  tolerant  of  these  uncertainties 
still  remains. 

C.  SUMMARY  OF  THIS  WORK 

In  this  dissertation  we  propose  a  method  for  reducing  the  number  of  real  para¬ 
metric  uncertainties  with  which  the  Hoo  control  algorithm  must  contend.  The  tech¬ 
nique  involves  grouping  uncert<un  parameters  which  affect  the  same  open  loop  eigen¬ 
value,  then  scaling  one  or  more  of  these  to  accomodate  the  effect  on  the  eigenvalue  of 
all  the  parameter  variations.  It  is  then  shown  that  the  resulting  controller  is  stable 
in  the  presence  of  all  the  uncertainties  while  still  meeting  the  performance  criteria. 

The  remainder  of  this  dissertation  is  organized  as  follows.  Chapter  II  provides 
a  summary  of  robust  multivariable  control,  including  the  use  of  the  structured  sin¬ 
gular  value  for  stability  analysis  and  the  state-space  equations  for  designing  an  Hoo 
controller.  Chapter  III  introduces  sensitivity  theory,  where  we  show  how  parameter 
uncertainty  affects  the  eigenvalues  of  the  open  loop  system.  The  technique  for  group¬ 
ing  parameters  which  affect  system  eigenvalues  similarly  and  scaling  one  or  more  of 
the  parameter  uncertainties  to  account  for  the  changes  caused  by  all  the  variations 
is  demonstrated.  This  technique  is  then  applied  to  the  design  of  control  systems  for 
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three  systems  in  Chapter  IV,  a  simple  four  state  system  and  two  missile  systems.  In 
Chapter  V  we  provide  conclusions  and  suggested  areas  for  further  research. 
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II.  ROBUST  MULTIVARIABLE  CONTROL 


Robust  multivariable  control  involves  the  design  of  feedback  controllers  for  sys¬ 
tems  with  specified  uncertainties  which  meet  established  performance  and  stability 
criteria.  For  purposes  of  this  work,  we  define  stability  and  performance  criteria  for 
systems  with  feedback  control  as  follows  (Dailey,  1990): 

•  Nominal  Stability:  stability  of  the  nominal,  unperturbed  system. 

•  Robust  Stability:  stability  of  the  system  for  a  given  set  of  perturbations  to 
the  nominal  model. 

•  Nominal  Performance:  the  nominal  system  meets  one  or  more  requirements 
for  closed  loop  performance. 

•  Robust  Performance:  the  system  meets  performance  requirements  and  re¬ 
mains  stable  for  a  given  set  of  perturbations  to  the  nominal  model. 

A  robust  Hoo  controller  is  a  feedback  controller  for  which  the  closed  loop  plant  exhibits 
robust  performance  as  defined  above.  The  concept  of  uncertainty  in  system  modeling 
is  presented  in  Section  A.  We  define  a  measure  of  robust  stability,  the  structured 
singular  value,  and  also  introduce  performance  robustness  in  Section  B.  Hoo  control 
theory  and  the  design  of  robust  Hoo  controllers  is  described  in  Sections  C  and  D. 

A.  UNCERTAINTY 

A  linear,  time  invariant  model  of  a  physical  plant  often  provides  an  adequate 
mathematical  representation  of  the  actual  system.  The  model  is  characterized  by 
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parameters  which  are  chosen  to  describe  the  system  at  a  given  operating  condition, 
such  as  when  the  system  is  at  rest  or  in  an  expected  steady-state  condition.  A  plant 
model  defined  in  this  way  is  termed  the  nominal  plant  model.  In  most  instances,  this 
nominal  model  is  but  one  of  an  infinite  number  of  models  required  to  actually  describe 
the  plant  at  any  moment  in  time.  Variations  in  the  plant  model  from  the  nominal 
conditions  are  m2aiifested  in  three  ways.  First,  a  number  of  pau’ameters  used  to 
describe  the  nomind  plant  may  be  uncertain,  or  known  only  within  some  measurable 
tolerance.  Examples  of  this  are  uncertain  pole  locations  or  damping  factors.  Secondly, 
certain  plant  dynamics  may  have  been  ignored  in  the  nominal  design,  such  as  high 
frequency  dynamics.  Lastly,  parameters  of  a  system  may  vary  as  operating  conditions 
of  the  physical  plant  change.  This  is  quite  common  in  aircraft  and  missile  control 
system  design,  where  changes  in  angle  of  attack  and  vehicle  weight  can  dramatically 
affect  the  aerodynamic  coefficients,  which  in  turn  affect  the  dynamics  of  the  plant. 
Hereinafter,  any  of  the  above  manifestations  which  cause  the  true  plant  to  vary  from 
the  nominal  plant  will  be  deemed  uncertainties. 

1.  Unstructured  Uncertainty 

Consider  the  plant  model  shown  in  Figure  2.1.  The  actual  plant  varies  from 
the  nominal  plant  by  the  uncertainty  Ai(s): 

G{s)  =  Go{s)[I  +  Ai{s)] 

where  Aj(s)  represents  a  frequency  dependent  uncertainty  in  the  plant  dynamics.  An 
uncertainty  represented  in  this  manner  is  termed  an  input  multiplicative  uncertainty 
and  it  can  used  by  a  designer  to  represent  unmodeled  dynamics  in  the  plant.  To 
this  uncertainty  we  associate  a  weighting  function  which  indicates  a  bound  on  the 
magnitude  of  the  uncertainty  as  a  function  of  frequency.  For  example  the  input  un¬ 
certainty  weighting  function  shown  in  Figure  2.2  indicates  a  50%  modeling  error  in 
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Figure  2.1:  System  with  input  multiplicative  uncertainty 

the  plant  below  10  rad/sec  and  a  100%  error  at  173  rad/sec  (Balas,  Doyle,  Glover, 
Packard,  Smith  1991). 

For  purposes  of  our  analysis  it  is  convenient  to  parameterize  the  uncertainty 
in  a  normalized  fashion.  As  a  way  of  illustration,  Figure  2.3  shows  how  the  uncer¬ 
tainty  Ai(s)  is  decomposed  into  a  frequency  dependent  weighting  function 
indicating  the  largest  singular  value  at  each  frequency  of  the  uncertainty,  and  a  nor¬ 
malized  uncertainty  A(s)  having  an  infinity  norm  not  exceeding  one.  The  system 
represented  in  this  way  is  said  to  be  in  standard  form,  with  the  plant  described  by  the 
nominal  plant  with  associated  weighting  function  designated  as  P(s).  The  transfer 
function  matrix  P(s)  now  has  two  inputs,  an  uncertainty  input  wj  and  an  exogenous 
input  w,  and  two  outputs,  an  uncertainty  output  j/j  and  a  reference  output  y.  The 
perturbation  block  A  can  vary  in  both  magnitude  and  phase,  constrained  only  by 
||A||oo  <  1.  Such  an  uncertainty  is  termed  an  unstructured  uncertainty. 

Another  type  of  uncertainty  is  the  inpu*  feedback  uncertainty,  depicted  in 
Figure  2.4.  This  type  of  uncertainty  can  be  used  by  the  designer  to  provide  gain  and 
phase  margins  in  design  of  the  Hoo  controller  as  described  in  Appendix  C.  Both  input 


Figure  2.3:  Uncertain  system  in  standard  form 
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Ai  0 

0  Aj 


Figure  2.5:  Standard  form  for  two  uncertainties 

multiplicative  and  feedback  uncertainties  will  appear  in  the  examples  of  uncertainty 
reduction  in  Chapter  IV. 

2.  Structured  Uncertainty 

Figure  2.5  shows  a  feedback  uncertainty  in  standard  form  combined  with 
an  input  multiplicative  uncertainty.  The  perturbation  block  A  now  has  two  blocks 
along  the  diagonal,  Ai  and  A2,  with  scaling  provided  such  that  ||Aj|oo  <  !•  In  this 
structure,  the  perturbation  block  is  no  longer  free  to  vary  in  all  elements,  as  the 
off  diagonal  blocks  are  now  zero.  Two  or  more  uncertainties  presented  in  this  way 
become  structured  uncertainties.  An  unstructured  uncertainty  is  a  special  case  of 
structured  uncertainty,  where  the  perturbation  consists  of  only  one  block. 

The  uncertainties  presented  above  describe  very  general  variations  in  the 
system  transfer  function.  Very  often,  however,  the  designer  has  some  knowledge 
about  the  specific  variations  of  the  parameters  in  the  plant  model.  The  toleramce 
to  which  the  parameters  have  been  mearsured  may  be  known  or  the  parameters 
may  vary  in  some  prescribed  manner  according  to  the  operating  conditions.  These 
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Figure  2.6:  System  with  parametric  uncertainty 

uncertainties,  characterized  by  real  scalar  variations  of  certain  parameters,  are  termed 
parametric  uncertainties. 

As  an  example.  Figure  2.6  shows  a  two  state  system  with  uncertainties  in 
the  pole  locations  plo  and  p2o.  The  uncertainties  and  62  are  defined  to  lie  in 
some  repon  constrained  by  61  €  (-ci,+ci]  and  62  €  [-62, +£2].  The  same  system 
is  shown  in  its  standard  form  in  Figure  2.7.  Although  the  uncertainties  and  62 
are  real  scalars,  the  perturbation  block  A  is  restricted  only  by  its  diagonal  structure 
and  infinity  norm;  Ai  and  A2  may  take  on  complex  values.  The  design  and  analysis 
techniques  presented  in  this  chapter  do  not  distinguish  between  real  and  complex 
perturbations,  therefore  real  parameter  uncertainties  are  treated  as  being  complex 
valued.  It  is  this  assumption  of  complex  blocks  in  the  perturbation  matrix  which  will 
subsequently  allow  us  to  reduce  the  number  of  parametric  uncertainties  needed  for  a 
robust  controller  design. 
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B.  ROBUST  STABILITY 


A  nominally  stable  feedback  system  can  become  unstable  in  the  presence  of 
uncertainties  in  the  nominal  plant  model.  When  designing  feedback  systems,  we 
need  a  stability  measure  other  than  just  stability  of  the  nominal  system,  and  therefore 
introduce  the  concept  of  robust  stability  as  a  measure  of  a  system’s  stability  in  the 
presence  of  uncertainties  in  the  plant.  In  this  section  we  will  define  a  measure  of 
robust  stability  and  describe  the  method  used  for  its  determination. 

1.  Stability  Robustness 

Consider  the  system  in  Figure  2.8  showing  the  uncertiiinty  block  A.  The 
plant  P(s)  contains  the  nominal  plant  and  uncertainty  weightings  such  that  ||  A(s)||oo  < 
1.  Associated  with  the  system  is  a  feedback  controller  K{s)  designed  on  the  basis  of 
the  nominal  plant  to  meet  given  specifications.  The  inputs  to  P{»),  wj,  w  and  u  are 
the  uncertainty  inputs,  exogenous  inputs,  and  control  inputs,  respectively.  The  out¬ 
puts  yj,  y  and  m  are  the  uncertainty  outputs,  reference  outputs  and  measurements, 
respectively. 

By  combining  P{s)  and  A'(s)  into  a  single  block  M{s),  the  system  can  be 
represented  as  in  Figure  2.9,  with  its  transfer  function  A/(s)  partitioned  as 

yd  _  A/ii(s)  Mi2{s)  1  [  1 

y  \  M2i{s)  M22(s)  J  [  u;  J 

Since  the  compensator  for  the  system,  A’(s),  is  designed  for  the  nominal  plant,  the 
system  Af(s)  is  nominally  stable.  To  guarantee  that  the  closed  loop  system  will 
remain  stable,  we  need  only  show  that  the  loop  contuning  the  perturbation  matrix 
is  stable  for  all  perturbations.  Due  to  the  normalization  of  the  perturbation  matrix  A, 
only  the  transfer  function  from  wj  to  yd,  namely  Mii(s),  need  be  analyzed.  Using  the 
stability  analysis  results  from  Appendix  B,  a  sufficient  condition  for  robust  stability 
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Figure  2.8:  Standard  form  for  an  uncertain  system  with  feedback  control 


Figure  2.9:  Model  for  robust  analysis 
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for  all  A  €  A  such  that  |iA||oo  <  1  is 


||Mn||oc<l.  (2.1) 

We  observe,  however,  that  the  only  constraint  this  result  places  on  the  perturbation 
matrix  is  a  bound  on  the  infinity  norm;  off  diagonal  elements  are  free  to  assume  a 
non-zero  vdue.  As  presented  earlier,  though,  the  perturbation  matrix  for  structured 
uncertainties  is  actually  block  diagonal.  The  stability  criterion  in  Equation  2.1  is  thus 
too  conservative  (Doyle,  1982).  In  the  next  section,  we  will  exploit  the  structure  of 
the  perturbation  matrix  and  develop  a  less  conservative  criterion  for  robust  stability. 

2.  The  Structured  Singular  Value 

An  alternative  to  Eiquation  2.1  for  determining  robust  stability  can  be  de¬ 
duced  from  Figure  2.9.  Stability  in  the  loop  containing  the  perturbation  requires 
that  for  all  frequencies 

det{/  -I-  MA)  ^  0  for  all  A  €  A  (2.2) 

where  A  is  the  set  of  all  possible  perturbation  matrices.  Doyle  has  defined  a  function 
which  can  be  used  to  define  the  necessary  and  sufficient  conditions  for  stability,  the 
Structured  Singular  Value  (SSV),  p(M)  (Doyle,  1982).  When  Equation  2.2  holds, 
this  function  has  the  property 

a{A)n{M)  <  1  (2.3) 

where  ^(A)  is  the  maximum  singular  value  of  the  perturbation  matrix  A.  The  SSV 
is  defined  as 

(M\  —  /  ®  if  no  A  €  A  solves  det(/  -f  MA)  =  0 

^  \  {min^g^(a(A)  |  det(/ -J- MA)  =  0))~*  otherwise 
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It  can  be  shown  that  the  SSV  has  the  following  properties  (Doyle,  1982): 

1.  n{aM)  =1  Q  I  ti{M) 

2.  p(/)  =  1 

3.  ti{AB)  <  a{A)fi{B)  A,B  complex  matrices 

4.  fi{A)  =  a(A) 

5.  DAD-^  =  A  'i  Dev 

6.  fi{DMD~^)  =  /x(A/)  D  a  nonsingular  diagonal  matrix 

7.  max  p{UM)  <  fi(M)  <  inf d^v  ^{DM D~^ ) 

where  p{X)  is  the  spectral  radius,  or  largest  eigenvalue,  of  X  and  the  two  sets  V  and 
Li  are  defined  as 

V  =  diag(di/ik,,</24j,---,d„/fc„),  d,  €  72; 

and  Ik  is  the  k  x  k  identity  matrix.  Unfortunately,  the  SSV  cannot  be  analytically 
determined  from  its  definition.  Numerical  methods  for  computing  the  SSV  exist, 
however  they  are  extremely  cumbersome  and  not  appropriate  for  most  design  work 
(de  Gaston  and  Safonov,  1988).  We  shall  therefore  use  Property  7  and  calculate  the 
bounds  for  the  SSV.  These  bounds  are  used  in  design  and  analysis  in  place  of  the 
actual  SSV. 


a.  Calculating  bounds  for  the  SSV 

Property  7  above  provides  upper  and  lower  bounds  for  the  SSV.  Doyle 
(Doyle,  1982)  has  proven  that  the  lower  bound  of  p,  given  by 

mjw  piMU)  <  p{M)  (2.4) 
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is  actually  an  equality.  This  problem  is  not  convex,  however,  and  local  maxima  can 
e  ast.  The  upper  bound  for  n  is  given  as 

ti{M)  <  mfj  a{DMD-' ).  (2.5) 

Determining  this  bound  is  a  convex  optimization  problem  and  it  has  one  minimum 
which  has  to  be  the  global  minimum.  The  minimum  value  of  a{DMD~^ )  can  be 
efficiently  computed  using  one  of  several  optimization  algorithms.  This  upper  bound 
has  been  found  to  be  generally  tight,  usually  within  5%  to  15%  of  the  SSV.  In  the 
remainder  of  this  dissertation,  we  will  use  the  SSV,  or  “mu"  to  describe  the  bound 
defined  by  Equation  2.5.  The  actual  value  of  ),  when  needed,  will  hereinafter  be 
referred  to  as  the  real  /z. 

3.  Performance  Robustness 

A  number  of  performance  criteria  can  be  defined  by  bounding  ||ry,i,lloo» 
where  Ty^u  is  the  transfer  function  from  the  set  of  exogenous  inputs  w  to  the  set  of 
error  measurements  y.  By  appropriately  weighting  the  output  vector  y,  we  define  the 
performance  criterion  such  that  it  is  bounded  by 

l|r„,||oo  <  1  (2.6) 

when  the  desired  performance  requirements  are  met  (Bibel  and  Malyevac,  1992). 
A  system  is  said  to  possess  nominal  performance  if  the  bound  in  Equation  2.6  is 
satisfied  for  the  nominal  plant.  The  system  is  said  to  possess  robust  performance 
if  the  system  remains  stable  and  meets  the  performance  criteria  in  the  presence  of 
all  uncertzunties.  Robust  performance  requirements  can  be  included  in  SSV  analysis 
by  adding  a  block  as  a  diagonal  term  of  the  uncertainty  matrix  A,  as  shown  in 
Figure  2.10.  The  reason  why  we  take  this  approach  is  because,  in  this  way,  robust 
performance  and  robust  stability  can  be  addressed  within  the  same  framework.  The 
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Figure  2.10:  Model  for  robust  performance 


performance  block  Ap  €  C'*’,  where  r  is  the  number  of  exogenous  inputs  and  q  is 
the  number  of  error  measurements.  The  augmented  uncertainty  inputs  and  outputs 
for  the  plant  are  now  defined  by 


Expressed  in  this  form,  the  test  for  robust  performauice  now  involves  computing 
the  SSV  for  the  entire  system  M{s),  rather  than  just  A/ii(s)  as  in  the  case  of  ro¬ 
bust  stability.  The  system  will  therefore  exhibit  robust  performance  provided  that 
/t(A/)  <  1.  It  is  worth  noting  that  performance  robustness  is  a  more  stringent  test 
than  stability  robustness. 

C.  Hoo  CONTROL  THEORY 

The  use  of  Hoo  control  theory  to  design  a  feedback  controller  for  a  multi-input 
multi-output  plant  is  introduced  in  this  section.  Hat,  design  involves  minimizing 
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Figure  2.11:  Model  for  //oo  controller  design 

the  peak  value  of  the  frequency  response  magnitude  of  selected  closed-loop  transfer 
functions.  The  efficacy  of  the  resulting  design  in  the  presence  of  uncertainties  and 
performance  criteria  is  then  analyzed  using  the  techniques  described  in  the  preceding 
sections. 

1.  Hoo  Optimization 

Consider  the  system  shown  in  Figure  2.11.  The  variables  are  defined  as 

follows: 

•  tn  €  'K-’"*  is  the  input  disturbance  vector,  including  uncertainty  inputs,  reference 
inputs  and  disturbances, 

•  ti  €  Ti”’*  is  the  control  vector  generated  by  the  controller, 

•  y  €  TV*  is  the  error  vector  containing  those  signals  we  desire  to  minimize, 

•  m  €  is  the  measurement  vector,  and  input  to  the  controller. 


19 


Since  the  system  is  assumed  to  be  linear  time  invariant,  by  partitioning  the  plant 
P(s)  we  can  write  the  outputs  in  terms  of  the  inputs  as 

Pu(s)  Pl2(s)  1  f  U;  1 

P21(S)  P«{S)  J  [  U  J 

or,  equivalently 

Jf  =  PjjW  +  P12U 

m  =  Pjiw  +  P32U 

Under  closed  loop  conditions  where  u  =  Km,  we  can  relate  the  disturbance  w  to  the 
error  signal  y  as 

y  =  [Pn  +  Pi2K{J  -  P22K)-^P2i)]w  =  P,{P,K)w. 

where  Pi{P,K)  is  termed  the  lower  linear  fractional  transformation  of  (P,  A').  Using 
the  above  definitions  for  w  and  y,  it  is  possible  to  put  a  number  of  design  problems 
into  the  form 

minimize  ||^/(P,A')||oo  (2.7) 

This  is  the  /foo-optimization  problem. 

2.  Solution  of  the  H^o  optimization  problem 

A  state-space  solution  to  Equation  2.7  is  summarized  here  (Doyle,  Glover, 
Khargonekar  and  Francis,  1988).  The  system  in  Figure  2.11  can  be  described  in  state 
space  form  as 

x{t)  =  Ax{i)  -1-  Biw{t)  -f-  B2u{t), 
y{t)  =  Duw{t)-\- Di2u{t), 

m{t)  =  C’2X(1)  -b  D2lXv{t)  -f-  D22U{t). 
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All  uncertainty  and  performance  function  weights  are  assumed  to  be  included  in  this 
plant.  The  plant  P{s)  now  has  a  trzmsfer  function  which  can  be  expressed  in  matrix 
form  as 


P(s)  = 


Di\  Di2 

Di\  D22 


+ 


Cl 

C2 


{sl-A)-'[Bi 


B2]. 


For  convenience,  the  plant  is  often  symbolically  expressed  as 


■  A 

Bi 

B2] 

P(s)  = 

Cl 

C2 

Dll 

D21 

Du 
D22  . 

which  is  not  a  transfer  function  matrix,  but  rather  a  compact  way  of  expressing  F(s). 
The  minimization  called  for  in  Equation  2.7  has  no  convenient  closed  form  solution. 


so  we  will  instead  solve  a  related  problem  and  find  a  controller  such  that 


mP,K)\\oo<% 


(2.8) 


where  7  is  a  predetermined  parameter.  This  will  involve  solving  two  Riccati  equations 
requiring  the  following  constraints  on  the  system  P{s): 

1.  (A,  82,02)  is  stabilizable  and  detectable. 

2.  rcink  D12  =  mj;  rank  D21  =  P2- 

3.  DJ2D12  =  /;  D21D21  =  /;  (For  Du  and  D21  full  rank,  scaling  can  be  used  to 
make  this  constraint  true.) 


4.  Dll  —  0;  D22  —  0- 


5.  rank 

6.  rank 


A  —  jul  B2 

Cl  Du 

A  —  jijjJ  Bi 

C2  D21 


=  n  +  m2 

=  n  +  p2 
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Denoting  the  solution  to  the  algebr&ic  Riccati  equation  A^X  ■\-X A  —  XRX  +  (J  =  0 
in  the  Hamiltonian  matrix  form 


V  D-  ^  -H 

we  now  define  two  matrices  Xoo  and  Yoo  as  the  solutions  to 


Xoo  =  Ric 


A  -B2BI 

-CfC,  -A^ 


K»  =  Ric 


-  CIC2 

-BiBj  -A 


Defining  three  intermediate  terms 


(2.10) 


z  =  {i-r^YX)-^-, 

F  =  -BlX\ 

L  =  -YCl, 


the  state-space  form  of  the  controller  K{s) 


Xc{t)  =  AcX^ii)  +  Bcm{t) 
u(f)  =  CeXc{t) 


is  defined  by  the  matrices 


Ac  =  A  +  ^-^BiBJX-\-B2F  +  ZLC2; 

Be  =  -ZL] 

Cc  =  F. 


The  following  conditions  must  be  met  for  a  compensator  K  satisfying  Equation  2.8 
to  exist: 


22 


1.  The  solution  to  Equation  2.9,  Xoo,  is  positive  semidehnite  and  the  associated 
Hamiltonian  has  no  imaginary  eigenvalues. 

2.  The  solution  to  Ek}uation  2.10,  Ko,  is  positive  semidehnite  and  the  associated 
Hamiltonian  has  no  imaginary  eigenvalues. 

3.  The  spectral  radius  p(A’ooVoo)  <  7^- 

The  controller  that  minimizes  the  cost  function  in  Equation  2.8  is  found  by  iteratively 
reducing  7,  with  the  initial  value  selected  large  enough  to  meet  the  above  criteria. 
The  value  of  7  is  reduced  using  the  bisection  method  until  no  improvement  in  7  is 
acheived  (Doyle,  Glover,  Khargonekar  and  Francis,  1988).  A  method  for  finding  the 
minimum  value  of  7  without  iteration  has  been  suggested  by  Chen,  Saberi  and  Ly 
(Chen,  Saberi  and  Ly,  1992). 

D.  D-K  ITERATION 

The  design  of  an  Hoo  controller  using  the  procedure  of  the  previous  section 
does  not  imply  robust  stability  or  performance.  The  algorithm  generates  a  controller 
which  only  guarantees  stability  of  the  nominal  closed  loop  system  and  minimizes 
||.Fi(P,  iif)lIoo-  The  resulting  closed  loop  system  must  still  be  tested  for  performance 
robustness  as  described  in  Section  B.  If  the  resultant  structured  singular  value  is 
not  less  than  one,  the  system  must  be  modified  and  another  controller  designed.  We 
now  introduce  the  technique  for  designing  a  closed  loop  system  that  has  performance 
robustness.  This  technique  uses  /foo  control  design  on  an  augmented  plant  and 
iteratively  modifies  the  plant  with  frequency  dependent  weighting  functions. 
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Figure  2.12:  Model  for  D-K  iteration 
1.  Modifying  the  Plant  With  D-Scales 

Recalling  from  Equation  2.5  that  the  actual  measure  of  robustness  is  the 
SSV  and  not  the  maximum  singular  value,  the  problem  stated  in  Equation  2.7  can 
be  more  accurately  posed  as 

minimize  \\DTi{P,  K)D~^  ||oo  (2-11) 

where  D  =  D{a)  is  a  frequency  dependent,  nonsingular  diagonal  scaling  function 
matrix.  It  can  be  shown  that  in  this  way  we  attempt  to  minimize  the  SSV  rather 
than  the  maximum  singular  value  (Dailey,  1990). 

The  appropriate  D-scales  are  selected  in  the  same  manner  as  in  the  Struc¬ 
tured  Singular  Value  problem  and  appended  to  the  plant  as  shown  in  Figure  2.12.  A 
new  /Too  controller  is  then  designed  for  the  augmented  plant  and  tested  for  robust¬ 
ness.  This  process  of  controller  design,  robustness  testing  and  D-scale  augmentation 
is  called  D-K  iteration  and  continues  until  the  maximum  value  of  the  SSV  is  less  than 
one  or  no  further  reduction  is  possible.  If  the  iteration  procedure  does  not  succeed 
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in  designing  a  robust  controller,  the  uncertunty  weights  or  performance  criteria  in 
the  plant  F(s)  must  be  modified  and  the  above  design  process  repeated. 
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III.  PARAMETRIC  UNCERTAINTY 
REDUCTION 

The  concept  of  plant  uncertainty  has  been  introduced  in  Chapter  II.  As  men¬ 
tioned  in  Chapter  I,  current  synthesis  techniques  may  fail  to  design  a  robust 
controller  when  all  the  known  uncertainties  are  included  in  the  plant  model.  In  this 
chapter  we  present  a  technique  for  reducing  the  number  of  parametric  uncertzunties 
used  in  the  design  of  a  robust  H^o  controller.  This  reduction  is  accomplished  by  using 
sensitivity  theory  to  identify  groups  of  parameters  whose  variations  affect  the  same 
open  loop  plant  eigenvalues,  and  then  quantifying  their  effect  on  the  eigenvalues.  The 
variation  in  one  or  more  of  these  parameters  is  then  scaled  in  a  manner  which  causes 
the  same  effect  on  the  eigenvalue  when  acted  upon  by  all  the  uncertainties.  An  Hoo 
controller  is  designed  for  the  system  with  the  reduced  number  of  uncertainties  and 
tested  for  robustness  with  the  system  containing  the  original  uncertainties. 

Sensitivity  theory  and  the  procedure  for  quantifying  the  effects  of  parameter 
variation  on  the  eigenvalues  of  the  open  loop  system  is  described  in  Section  A.  A 
technique  for  reducing  the  number  of  uncertainties  used  for  designing  a  robust  Hqo 
controller  is  presented  in  Section  B. 

A.  SENSITIVITY  THEORY 

In  this  section  we  present  a  discussion  introduced  by  Frank  (Frank,  1978)  of 
the  effect  of  real  pairametric  uncertaunties  on  the  nominal  plant.  We  will  be  most 
concerned  with  the  effect  the  pairameter  variations  have  on  the  system  eigenvalues. 
The  approach  taken  will  be  based  on  sensitivity  theory  by  developing  an  appropriate 
sensitivity  function.  This  function,  not  to  be  confused  with  the  sensitivity  transfer 
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function  of  the  closed  loop  system,  relates  the  change  in  the  parameter  of  interest  to 
specific  ch2U‘acteristics  of  the  plant.  We  first  present  the  underlying  theory  in  general 
form  and  then  define  the  sensitivity  function  relating  changes  in  plant  par2uneters  to 
changes  in  the  open  loop  eigenvalues. 

1.  General  Discussion 

Let  the  parameters  which  describe  a  system  be  represented  by  a  vector 
a  =  {qi,Q2, The  paraireter  vector  is  uncertain  around  a  nominal  vector 
Oo,  which  can  be  written  as 

a  =  Qo  +  Aot.  (3.1) 

For  linear,  time  invariant  systems,  the  plant  and  input  matrices  are  a  function  of  q, 
e.g-» 

X  =  A{q)x  -f  B{a)u.  (3.2) 

We  now  introduce  a  general  function,  (  =  ^(o),  which  characterizes  some  behavior 
of  the  system,  for  example  a  performance  index,  set  of  eigenvalues  or  the  state  of  the 
system. 

Defining  the  nominal  system  function  by  =  ^(Qo),  the  absolute  sensitivity 
function  is  defined  by 

=  i  =  l,2,...,m  (3.3) 

with  t  =  1,2, ...  ,m  relating  the  change  in  (  to  a  change  in  each  parameter  a,. 

For  small  perturbations,  the  total  perturbation  of  ^  due  to  changes  in  a  is  computed 

as 

=  (3.4) 

1=1 

We  next  discuss  the  method  for  determining  the  Sensitivity  Function. 
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a.  Calculating  the  sensitivity  function 

The  plant  parameter  vector  o  contains  elements  which  can  be  functions 
of  time  varying  parameters  themselves,  such  as  mass,  angle  of  attack,  etc.  The 
element  q,  can  then  be  represented  as  a,  =  Qi{pi,P2^---iPr)i  with  the  nominal 
vector  denoted  by  po-  The  goal  now  is  to  derive  the  sensitivity  function  for  ^  based 
on  changes  in  p: 

s',  =  ^1,.,  J  =  1.2,...,r  (3.5) 


Applying  the  chain  rule,  we  can  write 

dj  dQ2 

dpj  dQ\  dpj  ^  da2  dpj 


d(  dQ„ 

4*  Q  ^  )  J  —  1,2, 

Opj 


Evaluating  this  at  =  a(po),  and  using  the  definition  of  S^,(oo)  in  Equation  3.3, 
Equation  3.5  now  becomes 

5|,=|;si,^u.  ;  =  '.2 . r  (3.7) 

and  the  total  parameter-induced  uncertainty  can  be  expressed  as 

=  (3-8) 

J=1  i=l  ^P} 

The  quantity  doif  dp  j  must  be  defined  for  each  choice  of  the  parameter  vectors  q  and 
p.  In  the  next  section.  Equation  3.8  is  applied  to  the  eigenvalue  sensitivity  problem. 

2.  Eig:nvalue  Sensitivity 

We  now  consider  the  case  where  i  is  the  eigenvalues  of  the  system  juid 
Q  =  aij  the  elements  of  the  pl2uit  matrix  in  Elquation  3.2.  We  define  the  eigenvector 
sensitivity  as  follows: 


Given  an  LTI  system  with  dynamics  defined  by  Equation  3.2  with  A  = 
A{aij),  A  €  and  a,j  the  elements  of  the  matrix  A  having  eigenvalues  A*,  the 
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sensitivity  function  relating  perturbations  in  the  eigenvalues  to  perturbations  in  the 
elements  a,j  is  given  by 

“  'da~j 

wbe^e  Co  represents  the  nominal  values  of  the  matrix  A. 

To  develop  a  method  of  computing  this  sensitivity  function,  we  first  define 
the  right  and  left  eigenvectors  Vi  and  Wi  associated  with  the  eigenvalue  A,  of  the 
system  described  in  Equation  3.2.  The  right  eigenvectors  (i  =  1,2, ...,n)  are 
defined  by  Avi  =  A,t;j  and  the  left  eigenvectors  (t  =  1, 2, . . . ,  n)  by  A^Wi  =  A,i/\. 
Letting  A  equal  the  diagonal  matrix  of  eigenvalues  A,  aind  defining  t;  =  [fj,  tj2,  •  ■  • ,  Vn] 
and  w  =  [u;i,  102, /dots,  tUn],  the  matrix  form  of  the  eigenveci-or  relations  become 

•^1 

A  = 

An 


Av  =  vA 
Aw  =  toA 


The  right  and  left  eigenvectors  are  normalized  such  that  w^v  =  v^w  —  1  where  I  is 
the  identity  matrix.  Using  the  above  equations,  we  can  now  write  A  ^ls 

A  =  w^  Av. 


It  can  be  shown  that  the  sensitivity  matrix  5^*  is  directly  related  to  the  left  and 
right  eigenvectors  as 


(3.10) 


We  now  describe  the  method  for  determining  the  term  daij/dp.  We  define 
a  perturbed  parameter  vector 

P*  -  bl  +  Api,P2  +  Ap2, .  .  .  ,Pr  +  Apr] 
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For  small  perturbations  around  the  nominal  vector,  dcij/dp  can  be  approximated  by 
Aaij/Ap.  The  tenns  Aa^  are  the  elements  of  the  matrix 


AA  =  >4(ao)  —  >l(a) 

where  a  =  a(p*). 

For  one  parameter  variation,  Ap  =  Api,  and  AaijjAp  =  Aa,j/Ap(.  The  chain  rule 
is  used  to  find  the  eigenvalue  sensitivity  dXiddpi 


dpi  ~  hh  . 


(3.11) 


Recognizing  dXk/daij  as  the  eigenvalue  sensitivity  function  from  Equation  3.10  and 
assuming  small  perturbations,  Equation  3.11  becomes 


Apt  "  txh 


/=l,2,...,r 


(3.12) 


This  equation  provides  a  numerical  method  for  determining  the  sensitivity  of  each 
eigenvalue  to  any  parameter  in  the  system.  In  the  next  section  we  will  investigate 
the  problem  of  changes  in  the  eigenvalues  and  present  the  procedure  for  dealing  with 
perturbations  not  infinitesimally  small.  A  two  state  example  will  be  presented  which 
uses  the  methods  described  in  this  chapter  to  calculate  eigenvalue  sensitivity. 


B.  UNCERTAINTY  REDUCTION 

The  theory  for  calculating  the  sensitivity  of  system  eigenvalues  to  small  vari¬ 
ations  in  plant  parameters  was  derived  in  the  previous  section.  Not  all  the  system 
parameters,  however,  affect  the  eigenvalues.  In  this  chapter  we  deal  only  with  those 
parameters  that  affect  the  eigenvalues  of  the  open  loop  system,  and  the  example  in 
this  section  will  illustrate  that  each  parameter  typically  affects  only  a  small  num¬ 
ber  of  eigenvalues.  We  will  show  that  by  grouping  parameters  that  affect  the  same 
eigenvalue,  a  single  parameter  can  be  scaled  and  used  in  the  controller  design  to 
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accommodate  all  the  uncertainties  of  the  grouped  parameters.  In  subsection  1  the 
effects  of  parameter  variations  on  the  eigenvalues  are  illustrated,  and  the  case  of 
large  perturbations  will  be  included.  In  subsection  2  an  ex2unple  is  presented  which 
incorporates  the  techniques  previously  developed. 

1.  Calculating  the  Eigenvalue  Sensitivity 

In  Section  A,  the  sensitivity  of  the  kth  eigenvadue  of  a  system,  A*,  to  vari¬ 
ations  in  a  single  parameter  pi  was  defined  in  Equation  3.12.  In  order  to  select  an 
appropriate  parameter  for  scaling,  we  first  describe  the  manner  in  which  parameter 
variations  affect  pole  locations,  showing  that  the  maximum  eigenvalue  change  will 
occur  for  some  combination  of  the  maximum  parameter  variations.  Procedures  for 
dealing  with  perturbations  large  enough  that  the  Sensitivity  Function  changes  sig¬ 
nificantly  over  the  parameter  variation,  which  is  often  the  case,  will  be  discussed 
next. 


a.  Behavior  of  eigenvalues 

In  this  section  we  consider  the  perturbation  of  a  single  real  eigenvalue 
or  a  pair  of  complex  eigenvalues.  First,  consider  one  eigenvedue  of  a  system  in  the 
complex  plane.  The  eigenvalue  is  determined  by  the  equation 

A  +  a  =  0 

where  a  =  o(p),  the  parameter  vector  of  the  system.  The  variation  in  the  pole 
location  is 

AA  =  -Ao  (3.13) 

where  Aa  =  a{po)  -  q(Po  +  Ap).  Figure  3.1  describes  a  one  state  system  with  two 
uncertain  parameters,  P  =  Po-h  AP  and  K  =  Ko-¥  AK.  The  characteristic  equation 
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Figure  3.1:  Uncertain  one  state  system 


of  this  system  is 


A  +  ?K  =  0. 


(3.14) 


From  Equation  3.14,  it  is  apparent  that  the  maximum  eigenvalue  change  will  occur 
for  some  combination  of  the  maximum  parameter  perturbations. 

In  the  case  of  complex  roots,  the  eigenvalues  are  the  solutions  to 


A^  +  -f  -y  =  0 


(3.15) 


where  ^  =  /8(p)  and  7  =  7(p).  The  relationship  between  the  parameter  variations  and 
AA  is  best  illustrated  by  an  example.  Consider  an  uncertain  system  with  complex 
eigenvalues  as  shown  in  Figure  3.2.  For  the  time  being,  the  eigenvalues  are  assumed 
to  remain  complex  over  the  full  range  of  parameter  variations.  (This  constraint  will 
be  removed  later.)  In  this  particular  case,  the  characteristic  equation  of  the  system 
is  given  by 

A^  +  P,A  +  Pa  =  0  (3.16) 


with  roots 


(3.17) 


The  variation  in  the  real  part  of  the  pole  is  due  entirely  to  variations  in  Pi.  The 
variation  in  the  imaginary  part  of  the  eigenvalue  is  due  to  variations  in  both  Pi  and 
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Pi 


Figure  .2:  Uncertain  two  state  system 

P2.  Note  that  variations  in  P2  only  affect  the  imaginary  part  of  the  eigenvalue.  From 
Ex]uation  3.17,  it  is  clear  that  as  with  a  purely  real  pole,  the  maximum  eigenvalue 
change  will  again  occur  for  some  combination  of  the  maximum  parameter  variations. 
It  will  also  be  shown  that  the  sensivity  to  change  depends  on  the  direction  in  which 
we  move  the  pole,  and  in  some  instances  the  sensitivity  remains  constant  throughout 
the  change. 

In  general,  most  physical  systems  will  be  found  to  have  parameter  vari¬ 
ations  large  enough  that  the  eigenvalue  sensitivity  will  vary  considerably  over  the 
entire  rsuige  of  values.  In  uncertainty  reduction  we  will  be  most  concerned  with  the 
maximum  change  of  the  eigenvalue,  which  occurs  at  the  extreme  of  the  parameter 
variations.  In  practice,  sensitivities  for  all  parameters  are  calculated  first  for  the 
nominal  case.  Sensitivities  for  the  maximum  perturbation  case  are  calculated  only 
for  a  select  number  of  parameters,  as  described  in  the  following  section,  where  we  also 
show  that  the  point  of  maximum  eigenvalue  sensitivity  will  be  of  the  most  interest. 
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P2 


Figure  3.3:  Two  state  system  with  three  uncertain  parameters 

2.  Uncertainty  Reduction 

We  will  now  detail  the  steps  necessary  to  determine  the  eigenvalue  sensitiv¬ 
ity,  group  parameters  affecting  the  same  eigenvalues  and  scjde  a  parameter  variation 
to  accommodate  the  maximum  change  in  pole  location.  To  illustrate,  let  us  con¬ 
sider  a  numerical  example.  The  uncertain  two  state  system  shown  in  Figure  3.3  is 
characterized  by  the  parameters  PI,  P2,  and  K\,  where 

PI  =  Plo  +  ^Pl 
P2  =  P2o-><-Bp2 
A'l  =  K\o-\-6k\ 

The  nominal  parameters  and  associated  uncertainties  are 

Plo  =  8 

P2o  =  2 
Klo  =  100 

6pi  €  {-0.4 -Pl^, -1-0.4  •  PI  J 
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TABLE  3.1:  EIGENVALUE  SENSITIVITY 


Eigenvalue  Sensivitivy 

Evaluated  at 

AA/Ap 

1  AA/Ap  I 

Plo 

-50±;28.5 

pi 

plo  "F  Aplmaz 

-50±j48.1 

plo  Aplnto* 

-50±jl5.6 

p2o 

■EZM 

4.1 

p2 

p2o  +  Ap2fnax 

BBH 

3.3 

p2o  -  AP2„ax 

If 

5.9 

klo 

-0.07q:j0.04 

kl 

k\o  4"  Aklniox 

-0.07q:j0.03 

klo  Aklniox 

-0.07  q:  j0.05 

6p2  €  [-0.4P2„,+0.4P2„] 

6ki  €  [-0.2  A'l„,+0.2  A'g 

The  pole  locations  for  the  nominal  system  are  shown  in  Figure  3.4,  along  with  the 
pole  motions  for  the  above  parameter  variations.  Note  that  the  variations  in  PI 
and  A’l  result  in  real  and  imaginary  perturbation  of  the  pole  locations,  while  the 
variation  in  P2  results  in  purely  imaginary  motion  of  the  poles.  (This  is  consistent 
with  Equation  3.17.) 

Using  the  techniques  described  in  Section  A,  the  eigenvalue  sensitivities  from 
Equation  3.12  for  nominal  and  maximum  perturbation  cases  are  listed  in  Table  3.1. 
It  is  easy  to  see  how  the  sensitivity  changes  as  the  parameters  vary. 

As  discussed  in  Chapter  2,  the  parametric  uncertainties  can  be  arranged  in  a 
block  matrix  as  shown  in  Figure  3.5.  The  perturbation  block  A  has  real  scalar  values 
on  the  diagonal  and  all  other  elements  equal  to  zero.  The  only  restriction  placed  on 
the  perturbation  block  by  the  control  algorithm,  however,  is  that  the  frequency 
response  have  a  bounded  infinity  norm.  The  elements  in  A,  although  real,  may 
in  practice  assume  complex  values  without  affecting  the  controller  design  process. 
Without  loss  of  generality,  then,  we  can  express  each  element  in  A  as  a  bounded 
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-10 

real 


-5 


Figure  3.4:  Nominal  and  perturbed  pole  locations 
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magnitude  with  arbitrary  phase.  As  an  example.  Figure  3.6  shows  the  nominal  pole 
location  along  with  variations  in  PI  with  magnitude  )  API  |max  &nd  phase  varying 
from  0  to  2z.  The  pole  motion  describes  a  circle  around  the  nominal  pole  location. 
(The  circle  is  not  perfect  since  the  sensitivity  changes  with  parameter  variations  as 
described  above.) 

We  now  consider  variations  in  all  three  parameters  simultaneously.  For  three 
perturbations,  there  will  be  eight  (2^)  possible  pole  locations  when  evaluated  at  the 
scalar  extremes.  Figure  3.7  shows  the  nominal  and  perturbed  pole  locations.  The 
majcimum  AA  occurs  for  the  pole  at  —11.8  ±j2.4.  The  parameters  at  this  location 
are  described  by: 


PI  =  PU  +  APU„ 
P2  =  P2o-AP2,nax 
K1  =  Kh  +  AKlmas 


An  Poo  controller  designed  for  these  perturbations  should  then  accommodate  all 
variations  within  a  radius  of  |  AA^ai  I-  In  order  to  reduce  the  number  of  uncertainties 
used  for  Poo  design,  we  observe  that  the  variation  in  a  single  parameter  may  be 
increased  until  the  magnitude  of  the  eigenvalue  change  equals  that  of  the  extreme 
multiple  parameter  case.  In  the  above  example,  |  A  A  |mor=  10.99.  Table  3.1  suggests 
that  this  change  may  be  achieved  with  a  116%  variation  in  Pi,  calculated  as  follows. 
The  largest  eigenvalue  sensitivity  is  69.3  at  Pig,  +  APl„ior-  Using  this  to  provide  an 
initial  value  for  APlre,: 


Ap  = 

1  AA  [max 

AA/Ap 

1  AA  |max 

AP1„,  = 

AA/APl 

API 

0.1586  _ 

PI 

0.14 

10.99 

69.3 


=  0.1586 


(3.18) 
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Figure  3.6:  Pole  locations  for  arbitrary  phase  locations 
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Figure  3.7:  Nominal  and  Perturbed  Pole  Locations 
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Figure  3.8  shows  the  pole  location  associated  with  this  variation.  Note  that  116% 
variation  is  greater  than  that  required,  as  the  complex  poles  have  become  real.  (This 
is  not  unexpected  since  the  pole  sensitivity  increases  as  API  increases  in  the  pos¬ 
itive  direction.)  We  iterate  by  taking  successively  smaller  values  of  API  and  find 
that  the  required  |  AA  |max  is  found  to  occur  for  a  variation  of  95.5%  as  shown  in 
Figure  3.9.  An  Hoo  controller  may  now  be  designed  using  a  variation  of  95.5%  in  Pi. 
Uncertainties  of  40%  in  Pi  and  P2,  and  20%  in  A'l  should  then  be  accommodated 
in  the  design. 

The  parametric  uncertainty  reduction  process  can  be  summaurized  as  follows: 

1.  Calculate  the  eigenvalue  sensitivity  for  each  of  the  uncertain  parameters  at  its 
nominal  value. 

2.  Calculate  the  maximum  eigenvalue  change  for  each  eigenvalue  affected  by  the 
uncertain  parameters. 

3.  Group  parameters  affecting  the  same  eigenvalue  and  select  one  or  more  for  scal¬ 
ing.  (The  examples  in  Chapter  IV  will  provide  some  insight  into  this  selection 
process.) 

4.  Calculate  the  eigenvalue  sensitivity  for  the  maximum  parameter  perturbations 
and  scale  the  appropriate  parameter  uncertainties  to  accommodate  the  maxi¬ 
mum  eigenvalue  change  determined  above. 

5.  Use  the  scaled  uncertainties  in  the  Hoo  design  process  in  place  of  the  original 
uncertainties. 

6.  If  the  resulting  closed  loop  system  is  robust  in  the  presence  of  the  scaled  uncer¬ 
tainties,  verify  the  design  by  testing  the  controller  with  the  system  containing 
the  original  unsealed  uncertainties. 
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Figure  3.8:  Pole  location  With  116%  variation  in  pi 
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Figure  3.9:  Pole  location  with  95.5%  variation  in  pi 
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In  Chapter  IV,  we  use  three  ex2unples  to  demonstrate  the  above  uncertainty  reduction 
technique  in  the  design  of  a  robust  H^o  controller. 
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IV.  APPLICATIONS 


A  technique  for  reducing  the  number  of  parametric  uncertainties  needed  for  the 
design  of  a  robust  Hoo  controller  was  presented  in  Chapter  III.  In  this  chapter,  that 
technique  will  be  applied  to  three  dynamic  LTI  systems.  Parametric  uncertainty 
reduction  will  first  be  used  to  design  a  robust  Hoa  controller  for  a  simple  four  state 
system  with  five  uncertainties.  Next,  the  procedure  will  be  applied  to  the  design  of 
pitch  plane  autopilots  for  two  different  missile  models. 

A.  GENERIC  FOUR  STATE  MODEL 
1.  Model  Description 

An  uncertain  four  state  system  is  depicted  in  Figure  4.1.  The  system  is 
similar  to  that  used  in  the  example  in  Chapter  III  with  the  addition  of  two  more 
states  and  two  uncertain  parameters.  The  model  has  a  single  control  input  and 
feeds  back  one  state  and  an  error  signal  to  the  controller.  The  model  in  standard 
form  appears  in  Figure  4.2.  The  plant  has  two  exogenous  inputs,  the  reference 
and  a  disturbance  term.  (The  disturbance  has  been  added  to  fulfill  the  constraint, 
mentioned  in  Chapter  II,  that  Dji  be  of  full  rank.)  The  disturbance  has  a  nomind 
weight  of  Wditi  =  0.001,  There  are  five  uncertainties  and  two  performance  weights. 
The  controller  generates  a  single  command  input.  All  five  parameters  in  the  plant 
are  assumed  to  be  uncertain  around  a  nominal  value.  Table  4.1  lists  the  nominal 
value  for  each  parameter  and  its  associated  uncertainty. 
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Figure  4.1:  Four  state  system 


TABLE  4.1:  NOMINAL  SYSTEM  PARAMETERS  AND  UNCERTAINTIES 


Parameter 


Uncertainty 


a.  Performance  Weighting  f\inctions 


We  define  two  performance  weighting  functions  for  the  system,  a  weight 
on  the  error  between  the  reference  and  the  system  output  and  a  weight  on  the  control 
input.  The  control  weight,  VK»n<,  is  set  arbitrarily  small  at  0.001.  We  include  the 
weighted  control  as  an  output  only  to  fulfill  the  constraint  that  the  matrix  term  D12 
be  of  full  rank  as  mentioned  in  Chapter  II.  The  error  weighting  function  Werr,  plotted 
in  Figure  4.3,  is  a  low  pass  filter  with  low  frequency  gain  of  50,  a  gain  crossover  of  2 
rad/sec  and  a  high  frequency  gain  of  0.125.  This  choice  of  weighting  function  should 
result  in  a  steady-state  error  of  less  than  2%,  closed  loop  time  constant  of  0.5  and 
minimal  overshoot  (Bibel  and  Malyevac,  1992).  The  transfer  function  of  the  error 
weight  is 


W„r  = 


50(0.025s  +  1) 


(4.1) 


(10s -1-1) 

The  goal  is  to  design  an  Hoo  controller  which  meets  stability  and  robustness  criteria 
taking  into  account  all  hve  uncertainties. 


2.  Hoc  Controller  Design 


We  first  attempt  to  design  a  controller  by  including  all  five  uncertainties  in 
the  plant  model.  The  Hoo  controller  design  procedure  of  Chapter  II  was  used  in  an 
attempt  to  generate  a  robust  controller.  The  mu  plot  in  Figure  4.4  shows  the  SSV 
for  robust  performance  and  stability  and  the  principle  gains  for  nominal  performance 
after  three  D-K  iterations.  It  is  clear  that  the  closed  loop  plant  does  not  exhibit 
nominad  or  robust  performance.  Further  iterations  yielded  no  improv  ement.  We  next 
apply  the  uncertainty  reduction  technique  to  the  system  and  again  attempt  to  design 
a  robust  controller,  this  time  with  a  reduced  number  of  uncertainties. 
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Figure  4.4:  Robust  analysis 
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3.  Uncertainty  Reduction 

The  eigenvalues  of  the  nominal  system  are 

>1.2  =  -7±jl2.3; 

>3  =  -32.0; 

>4  =  0. 

The  eigenvalue  sensitivity  for  all  the  plant  parameters  are  shown  in  Table  4.2.  We 
observe  that  >1,2  are  affected  only  by  parameters  Pl,/*2,  and  Kl,  A3  by  parameters 
P3  and  ^*2,  while  the  pole  at  the  origin  remdns  unaffected.  The  maximum  eigenvalue 
changes  for  all  combinations  of  parameter  variations  was  found  using  the  procedure 
in  Chapter  III  to  be 

|AA,.2Uax  =  10.99; 
lAAaUo*  =  21.76. 

We  select  uncertainties  P\  and  PZ  for  scaling  and  apply  Equation  3.18  to  obtain  an 
initial  value  for  the  scaled  uncertainties: 

API  =  1.16  ♦PU 
APS  =  0.68  ♦P3„ 

The  eigenvalue  sensitivity  for  P3  remains  constant  over  the  range  of  variation,  and 
iteration,  described  in  Chapter  III,  is  not  required.  The  eigenvalue  sensitivity  for 
parameter  PI  varies,  however,  and  we  iterate  to  obtain  the  actual  variation  in  Pi 
necessary  to  cause  the  required  eigenvalue  change, 

APlre,  =  0.955  ♦  Plo 
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TABLE  4.2:  EIGENVALUE  SENSITIVITY 


Eigenvalue  Sensivitivy  | 

Parameter 

Evaluated  at 

AXlAp 

1  aa/Ap  I 

Plo 

-50±j28.5 

57.3 

PI 

Pl^  ^Plm»* 

-50±;48.1 

69.3 

PU  -  APUa. 

-50±jl5.6 

52.4 

P2o 

■oiinii 

4.1 

P2 

P2c  +  AP2max 

3.3 

-7±jl2.3 

P2„  -  AP2„„, 

0tj5.9 

5.9 

KU 

K1 

Klo  •¥  AKlmaz 

KU  -  ^Kl„az 

P3 

PZo 

0 

K2 

K2o 

0 

0 

PI 

Plo 

0 

0 

P2 

P2o 

0 

0 

K1 

KU 

0 

0 

PZo 

4 

-32.0 

P3 

PZo  +  ^PZynaz 

-4-f;0 

4 

PZo  ~~  AP3»no* 

-4  +  jO 

4 

K2o 

-8-l-jO 

8 

K2 

K2o  +  AK2maz 

-8-l-;0 

8 

K2o  -  AK2„oz 

-8-|-;0 

8 

PI 

PU 

0 

0 

P2 

P2o 

0 

0 

0 

Kl 

KU 

0 

0 

PZ 

PZo 

0 

0 

K2 

K2o 

0 

_ 2 _ 
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4.  Controller  Design  for  Reduced  Uncertainties 


Using  only  the  scaled  uncertainties  in  P\  and  P3,  it  was  possible  to  design 
an  Hoa  controller  exhibiting  desired  robustness  properties  in  the  presence  of  the  scaled 
uncertainties.  The  mu  plot  in  Figure  4.5  displays  the  SSV  for  robust  performance 
and  stability,  and  the  principle  gains  for  nominal  stability.  All  plots  remain  below 
one,  indicating  a  robust  design.  Our  interest,  however,  is  in  the  original  plant  with 
five  uncertainties.  Using  the  controller  designed  for  the  scaled  uncertainty  plamt, 
robust  analysis  of  the  complete  plant  resulted  in  the  mu  plot  of  Figure  4.6.  It  is 
clear  that  the  system  meets  robust  stability  and  nominal  and  robust  performance 
specifications. 

We  next  examine  the  time  response  of  the  final  closed  loop  system  to  verify 
that  performance  objectives  are  met.  Figure  4.9,  showing  the  system  response  to  a 
unit  step  input,  reveals  a  rise  time  of  0.23  seconds,  a  maximum  overshoot  of  4.4%  and 
a  steady  state  error  of  1.73%.  This  confirms  that  the  performance  goals  of  minimal 
overshoot  and  less  than  2%  steady-state  error  are  satisfied. 

5.  Controller  Design  Ignoring  Some  Uncertainties 

It  is  interesting  to  investigate  the  effect  of  controller  design  ignoring  several 
of  the  uncertainties  in  the  system.  We  may  be  tempted  to  use  this  approach  when 
the  inclusion  of  all  uncertainties  results  in  a  failure  to  design  a  robust  controller,  as 
in  subsection  2.  Considering  only  the  40%  uncertainties  in  PI  and  P3,  and  ignoring 
uncertainties  in  P2,  K\  and  K2,  an  //„  controller  was  designed.  When  tested 
against  the  plant  with  all  five  uncertainties,  however,  the  SSV  exceeds  one  as  shown 
in  Figure  4.8,  implying  that  the  robust  performemce  objective  is  not  met.  It  is  clear 
that  simply  ignoring  uncertainties  may  often  result  in  a  closed  loop  plant  which  lacks 
robustness. 
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frequency  (rad/sec) 


Figure  4.5:  Robust  analysis  for  scaled  uncertainty  plant 
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mu 


Figure  4.6:  Robust  analysis  for  plant  with  all  uncertainties 
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Figure  4.7;  Step  response 
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mu 


Figure  4.8:  Robust  analysis  for  plant  ignoring  some  uncertainties 
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Figure  4.9:  Low-drag  ramjet  model 

B.  COMBUSTOR  CONTROLLED  MISSILE 

We  next  apply  the  uncertainty  reduction  technique  to  a  pitch  plane  autopilot 
for  a  low-drag  ramjet  missile  model.  This  model,  shown  in  Figure  4.9,  is  controlled 
by  a  moveable  combustor  section  and  has  no  external  control  surfaces  on  the  missile. 
The  model  has  five  uncertain  parameters  affecting  two  real  poles  and  a  set  of  complex 
poles.  The  linearized  model  of  this  missile  was  developed  by  the  Naval  Air  Warfare 

Center,  China  Lahe,  CA  (Robins,  1992). 

✓ 

1.  Missile  Model 

The  block  diagram  of  the  missile  is  shown  in  Figure  4.10.  The  input  to  the 
system  is  a  commanded  combustor  defiection  angle,  6c  in  radians.  The  outputs  of  the 
system  are  missile  acceleration  t],  and  pitch  rate  q.  The  missile  model  was  linearized 
at  one  operating  altitude  and  velocity.  The  varying  factors  affecting  the  plant  are 
angle  of  attack  and  missile  mass.  Seven  parameters  are  affected  by  these  changes, 
Jiay  Hsi  Za,  Zg,  Mat  Mg  and  Jg.  Of  these,  Ha  &nd  Hg  do  not  affect  the  eigenvalues. 
Table  4.3  lists  the  nominal  value  for  each  pareunete''  along  with  uncertainties  in  the 


TABLE  4.3:  NOMINAL  SYSTEM  PARAMETERS  AND  UNCERTAINTIES 


Parameter 

Uncertainty 

h.  =  201 

Mi,  =  341.59 

Zi,  =  0.2165 

Da  =  1.17 

Ha  =  11,595 

Ic  =  14.5 

Lafg  =  0.1739 

Bm  =  1523 

varying  parameters  listed  above. 

a.  Performance  Weighting  Functions 

For  this  system  we  define  two  performance  weighting  functions,  one 
weighting  the  error  between  the  reference  input  and  output  acceleration  and  one 
weighting  the  combustor  deflection  angle.  The  error  weight  is  a  low  pass  filter  shown 
in  Figure  4.11.  It  has  a  DC  gain  of  50,  corresponding  to  a  2%  steady-state  error. 
The  combustor  deflection  weight  is  set  equal  to  the  inverse  of  the  maximum  allowable 
deflection  of  10  degrees. 

2.  Hqo  Controller  Design 

An  attempt  was  made  to  design  a  robust  Hoo  controller  including  all  five 
uncertainties  in  the  plant.  The  SSV  plot  of  Figure  4.12  reveals  that  the  process  failed 
to  p'oduce  a  controller  exhibiting  robust  performance.  We  next  apply  the  parametric 
uncertainty  reduction  technique  and  again  attempt  to  design  a  robust  controller. 
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frequency  (rad/sec) 


Figure  4.12:  Robust  analysis 
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3.  Uncertainty  Reduction 


The  eigenvalues  of  the  nominal  system  are 

Xia  =  -36.76  ±j62.7 
A3  =  27.35 
X4  =  -27.27 
As  =  -257.11 

and  their  eigenvalue  sensitivities  are  listed  in  Table  4.4.  Eigenvalue  A5  is  not  affected 
by  any  of  the  varying  parameters  and  has  not  been  included.  Following  the  procedures 
in  Chapter  III,  the  maximum  eigenvalue  changes  for  all  possible  combinations  of 
parameter  variations  was  found  to  be 

I  AAi^2  Irnoi  =  1-63 
lAAaUa*  =  4.47 
lAA^Uax  =  5.10 

We  select  the  parameters  for  scaling  by  exeunining  Table  4.4.  Parameter  1/h,  al¬ 
though  affecting  all  four  eigenvalues  with  large  sensitivities,  is  sufficiently  small  that 
a  variation  of  almost  200%  would  be  required  to  accommodate  the  change  in  A3. 
Likewise,  parameters  Zq  and  Zs  have  relatively  small  nominal  values  and  sensitivi¬ 
ties.  We  therefore  choose  to  scale  Ms  to  satisfy  the  change  in  Ai,2  and  Mo  for  the 
change  in  A3  and  A4.  The  initial  scaling  values  are  found  as  follows: 


AMs  = 


I  AA^.2  I  max 

AAi,2/AM{ 


1.63 

0.006 


271.67 


AM,  271.67 
Ms  341.59 


63 


TABLE  4.4:  EIGENVALUE  SENSITIVITY 
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AAf» 


I  AA3  |fnttx 

AAa/AM^ 


4.47 

0.021 


212.9 


AMa  212.9 
Ma  ~  760.79  ~ 


AMa 


I  AA4  Itngg 

AXJAMa 


5.10 

0.021 


242.9 


AMa  242.9  ^ 

A/a  “  760.79  " 

By  iterating,  the  required  scaled  uncertainties  are  found  to  be  80%  in  Ms  and  34% 
in  A/a.  These  scaled  uncertainties  are  next  used  to  design  a  robust  Hoc,  controller. 

4.  Controller  Design  for  Reduced  Uncertainties 


Using  only  the  scaled  uncertainties  in  parameters  Ms  and  A/a,  a  robust  Hoo 
controller  has  been  designed.  The  plot  of  the  SSV  is  shown  in  Figure  4.13.  The  design 
exhibits  nominal  performance  and  robust  performance  and  stability  in  the  presence 
of  the  two  scaled  uncertainties.  The  controller  when  tested  against  the  plant  with  all 
five  uncertainties  yields  the  plot  in  Figure  4.14.  The  system  again  exhibits  nominal 
and  robust  performance  and  robust  stability. 

The  time  response  of  the  final  closed  loop  system  was  analyzed  to  verify 
that  performance  criteria  were  met.  The  step  response  of  the  system,  shown  in 
Figure  4.15,  shows  an  overdamped  system  with  a  rise  time  of  0.35  seconds  and  no 
overshoot.  Steady  state  error  was  determined  to  be  0.5%.  The  initial  transient  spike 
is  a  result  of  the  torque  caused  by  the  movement  of  the  combustor  section. 
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Figure  4.13;  Robust  analysis  for  scaled  uncertainty  plant 
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Figure  4.14:  Robust  analysis  for  plant  with  all  uncertainties 
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Figure  4.16:  Tail  fin  controlled  missile 

C.  TAIL  FIN  CONTROLLED  MISSILE 

The  parametric  uncertainty  reduction  technique  is  next  applied  to  the  pitch 
plane  autopilot  of  a  conventional  tail  fin  controlled  missile.  The  missile  is  depicted  in 
Figure  4.16.  This  model  has  three  uncertain  parameters  affecting  one  set  of  complex 
eigenvalues.  One  of  the  parameters  varies  270%,  however,  causing  the  eigenvalues  to 
shift  from  complex  to  real  over  the  range  of  uncertainty.  A  linearized  model  from 
the  Naval  Surface  Warfare  Center,  Dahlgren,  VA  was  used  for  analysis  (Bibel  and 
Stalford,  1991). 

1.  Missile  Model 

The  block  diagram  of  the  missile  is  shown  in  Figure  4.18.  The  input  to  the 
system  is  a  commanded  tail  fin  deQection  angle,  5e.  The  system  outputs  are  missile 
acceleration,  q,  and  pitch  rate,  9,  from  the  accelerometer  and  gyro,  respectively. 
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Flexible  Body  Dynamics 


Figure  4.17:  Conventional  missile  block  diagram 
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TABLE  4.6;  FLEXIBLE  BODY  MODE  PARAMETERS 


a.  Nominal  System  Parameters  and  Uncertainties 


The  linear  model  has  been  developed  for  one  pju’ticular  altitude  and 
velocity.  Factors  affecting  the  plant  parameters  are  angle  of  attack,  a,  and  pitch  rate, 
q  which  cause  changes  in  the  parameters  Mo,  Za  and  Mg.  The  nominal  parameters 
and  associated  uncertainties  are  listed  in  Table  4.5.  The  actuator  is  modeled  as  a 
second  order  transfer  function  with  a  damping  of  (act  =  0.7  and  a  natural  frequency 
of  uact  =  188.5  rad/sec.  The  actuator  has  a  fin  deflection  limit  of  40  degrees  and  a 
fin  rate  limit  of  300  degrees/sec.  The  input  to  the  actuator  is  the  command  input, 
6e.  The  model  includes  flexible  bending  modes  which  are  modeled  as  second  order 
transfer  functions.  The  first  three  bending  modes  have  been  included,  and  the  nth 
mode  effect  on  acceleration  and  pitch  rate  is  given  by 


9(s)  _ 

1 

‘K,.  ■ 

.■JW  J/6‘ 

S*  +  2(fbUf„3  +  ul 

The  data  for  the  flexible  modes  is  shown  in  Table  4.6.  The  controller  acts  on  an 
accelerometer  and  pitch  rate  gyro,  both  modeled  as  second  order  devices  with  a 
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damping  of  0.7  and  natural  frequencies  u>acc  =  377  rad/sec  and  ucyro  =  500 
rad/sec.  The  model  includes  an  input  feedback  uncertainty  with  a  constant  value 
of  0.2,  and  a  second  order  input  multiplicative  uncertainty  with  a  triinsfer  function 
given  by 

^  (l/100)(0.049s  +  1)^ 

(9.805e-‘s  +  l)* 

b.  Performance  Weighting  Functions 


For  this  plant  we  define  four  performance  weighting  functions.  Two  of 
these  are  weights  on  the  fin  deflection  and  fin  rate,  set  equal  to  the  inverse  of  the 
allowable  values  of  40  degrees  and  300  degrees/second,  respectively.  The  performance 
weight  on  the  error  between  the  reference  input  and  output  acceleration  is  a  first  order 
transfer  function  described  by 


I00(0.0258s  +  1) 
10.327s -f  1) 


(4.2) 


The  fourth  performance  function  is  a  weight  on  the  pitch  rate.  It  is  a  high  pass  trans¬ 
fer  function  designed  to  attenuate  the  effects  of  unmodeled  high  frequency  dynamics 
and  sensor  noise.  It  is  described  by 


^  (l/400)(0.025s  +  1)" 
’  (1.78e-'‘s-|- 1)2 


(4.3) 


The  goal  is  to  design  a  robust  Hx,  controller  taking  into  account  all  the  uncertainties 
and  performance  weights. 


2.  Hoo  Controller  Design 


An  attempt  was  made  to  design  a  robust  H^o  controller  for  the  pleint  with 
all  three  uncertainties.  Figure  4.18  shows  the  robust  analysis  results  after  four  D-K 
iterations.  The  closed  loop  system  does  not  possess  robust  or  nominad  performance. 
Uncertainty  reduction  techniques  are  next  applied  to  the  model  in  an  effort  to  obtain 
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mu 


Figure  4.18:  Robust  analysis 
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TABLE  4.7:  EIGENVALUE  SE^’SITIVITY 


Eigenvalue  Sensivitivy 

Eigenvalue 

Parameter 

Evaluated  at 

AA/Ap 

1  AA/Ap  I 

Ma^ 

0  ±  jO.073 

Ma 

Mae  "f”  ^Mamex 

0  ±  jO.038 

Mae  -  AA/a„„ 

±0.056  ±j0 

0.5  4:;0.02 

0.5 

-4.97±j7.5 

2a 

2ae  "h  ^2onex 

0.5Tj0.05 

0.5 

2ae  ~  ^2a„ex 

0.5t;0.01 

0.5 

M,e 

0.5  ±  j0.02 

WSSM 

M, 

M,^  -f  AA/,„„ 

0.5  ±  j0.03 

Mqo 

0.5±;0.07 

0.5 

a  robust  controller. 


3.  Uncertainty  Reduction 


The  eigenvalues  of  the  open  loop  system  are 


•^1,2  = 

-4.97±j7.5 

-^3.4  = 

-38.5±jl923 

II 

tD 

-7.36  ±  j368 

A?, 8  = 

-18.7±j937 

^9,10  = 

-132±jl35 

Only  Ai,2  are  affected  by  the  variations  in  parameters  A/q,  Zq  and  A/,.  The  eigenvalue 
sensitivities  for  these  parameters  is  listed  in  Table  4.7.  The  parameters  for  scaling  are 
selected  by  observing  that  the  large  variation  in  A/q  causes  the  eigenvalues  to  vary 
from  complex  to  real.  While  complex,  the  largest  eigenvalue  change  is  |  AAi,2  1=  7.5, 
the  distance  from  the  nominal  location  to  the  real  axis  This  change  is  already  acco¬ 
modated  by  the  first  100%  variation  in  Ma-  On  the  real  axis  with  all  three  parameters 
varying,  the  eigenvalues  have  extreme  values  of  4-5.10  and  -15.1.  The  variation  in  Mq 
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alone  causes  the  eigenvalues  to  move  to  -1-4.03  and  -13.96.  The  ^u^ditional  eigenvaJue 
change  caused  by  variations  in  Za  and  can  be  accommodated  by  increasing  the 
variation  in  Mq  to  300%  and  leaving  unchanged  the  15%  vcu-iation  in  Zq.  An 
controller  is  next  designed  for  a  system  with  these  two  uncertainties  used  in  place  of 
the  uncertainty  in  Mq. 

4.  Controller  Design  for  Reduced  Uncertainties 

Using  only  the  scaled  uncertainties  in  Ma  and  Z^,  a  robust  Hoo  controller 
was  designed.  The  SSV  plot  for  the  closed  loop  system  is  shown  in  Figure  4.19.  When 
tested  against  the  system  incorporating  all  three  uncertainties  at  their  actual  values, 
the  plot  of  Figure  4.20  results.  The  system  meets  the  robustness  specifications. 

The  time  response  of  the  final  closed  loop  system  was  analyzed  to  ZLSsess 
performance.  The  step  response  of  the  system  is  show  in  Figure  4.21.  The  response 
shows  an  overdamped  system  with  a  rise  time  of  approximately  0.15  seconds  and 
no  overshoot.  The  error  at  1  second  is  1.4%,  and  has  decreased  to  near  the  1% 
specification  by  2  seconds. 
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Figure  4.19;  Robust  analysis  for  scaled  uncertainty  plant 
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frequency  (rad/sec) 


Figure  4.20:  Robust  audysis  for  plant  with  all  uncertainties 
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V.  SUMMARY  AND  CONCLUSIONS 


A.  PROBLEM  STATEMENT 

Current  Hoo  synthesis  techniques  are  often  unable  to  provide  a  robust  controller 
design  when  all  parametric  uncertainties  in  a  system  are  included  in  the  plant  model. 
The  plant,  in  effect,  becomes  overconstrained  and  the  design  algorithms  fail.  This 
does  not,  however,  preclude  the  existence  of  a  robust  Hoo  controller,  rather  it  indicates 
current  methods  are  not  effective  in  many  practical  design  problems.  This  work 
has  focused  on  a  technique  which  allows  the  designer  to  use  existing  Hoo  synthesis 
theory  and  commercially  available  software  to  produce  a  controller  which  is  robust 
in  the  presence  of  all  quantified  p<irametric  uncertainties  in  the  plemt.  The  procedure 
involves  grouping  uncertain  parameters  which  affect  the  same  open-loop  eigenvalue. 
The  Vciriation  of  one  or  more  of  the  designated  parameters  is  scaled  to  provide  for 
the  same  eigenvalue  change  as  that  caused  by  all  the  parameters  at  the  extreme 
values  of  their  uncertainties.  An  Hoo  controller  which  is  robust  in  the  presence  of  the 
scaled  uncertainties  is  then  shown  to  be  robust  in  the  presence  of  the  actual  plant 
uncertainties. 

B.  CONTRIBUTIONS 

As  a  result  of  this  research,  two  original  contributions  were  made  that  pertain 
to  robust  controller  design: 

•  Development  of  a  technique  which  reduces  the  number  of  parametric  uncer¬ 
tainties  required  for  the  design  of  a  robust  Hoo  controller  that  accommodates 
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all  parametric  uncertainties  in  the  plant,  when  design  attempts  using  all  un¬ 
certainties  have  failed.  This  technique  may  also  be  used  to  reduce  the  number 
of  uncertainties  used  even  when  a  robust  controller  can  be  found  for  all  uncer¬ 
tainties,  thus  reducing  computation  time. 

•  Procedures  for  applying  the  technique  to  missile  autopilots,  providing  for  the 
design  of  a  robust  controller  taking  into  account  not  only  model  uncertainties 
but  also  parameter  variations  caused  by  changes  in  flight  conditions.  These 
variations,  often  quite  extreme,  have  until  now  required  conservative  autopilot 
designs  which  sacrificed  performance  to  maintain  stability. 

In  short,  this  research  has  provided  designers  with  a  new  tool  for  finding  a  robust 
Hoo  solution  in  specific  cases  when  current  design  methods  fail.  Development  of  this 
technique  has  revealed  several  other  areas  requiring  further  study.  These  include: 

•  Complex  scaling  of  a  parametric  variation  rather  than  real  scaling  as  described 
herein  may  allow  for  the  use  of  a  smaller  scale  factor  in  robust  design,  and 
therefore  increase  the  number  of  problems  for  which  a  solution  can  be  obtained. 

•  Performance  functions  in  missile  autopilot  models  are  often  the  dominimt  con¬ 
straints  in  controller  design.  A  technique  for  reducing  the  number  of  perfor¬ 
mance  constraints  similar  to  the  technique  for  parametric  uncertainty  reduction 
may  also  increase  the  number  of  problems  for  which  a  solution  can  be  obtained. 
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APPENDIX  A 

Transfer  Function  Norms 

Structured  Singular  Value  (SSV)  analysis  makes  extensive  use  of  transfer  func¬ 
tion  norms  in  defining  the  underlying  theory  for  stability  analysis.  We  present  here 
a  brief  discussion  of  the  2-norm  and  oo-norm  along  with  some  useful  properties  (Ma- 
ciejowski,  1989). 

Given  a  proper  transfer  function  matrix  G(s)  with  no  poles  on  the  imaginary 
axis,  the  2-norm  of  G(s)  is  defined  as 

IIGIIa  £ 

The  infinity  norm  is  defined  as 

llGlloo  =  supa(G(iu;)) 

W 

with  a  denoting  the  largest  singular  value.  These  norms  satisfy  the  following  prop¬ 
erties: 

1-  ll^ll  ^  0  with  ((Gil  s=  0  if  and  only  if  G  =  0. 

2.  liaGII  =  |a|||G||  Va€C. 

3.  ||G-f^|l<l|G||  +  ||/f||. 

The  oo-norm  also  satisfies 


llG/flU  <  IlGlUlli/lloo 


This  inequality  is  not  satisfied  by  the  2-norm. 
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APPENDIX  B 

Robust  Stability  Analysis 


We  present  in  this  Appendix  a  discussion  of  robust  stability  analysis  using  the 
small  gain  theorem  (Burl,  1993).  An  uncertain  system  with  a  feedback  controller 
designed  for  the  nominsJ  plant  can  be  depicted  as  in  Figure  B.l.  By  combining  the 
nominal  plant  and  controller  into  one  block  Af(s),  we  can  redraw  the  system  as  in 
Figure  B.2.  Writing  the  outputs  of  A/(s)  in  terms  of  the  inputs: 

j/d  ]  _  [  A/„(s)  A/i2(s)  ]  f  Wrf  1 
y  \  Afjijs)  Afj2(s)  J  [  UJ  J 

Figure  B.3  shows  the  system  expanded  to  show  the  subsystems  of  Af(s),  The  nominal 
system  M(s)  is  stable  since  K{s)  is  designed  for  the  nominal  plant.  The  entire  system 
will  therefore  be  internally  stable  provided  that  the  loop  containing  A(s)  is  stable, 
since  this  is  the  only  possible  source  of  instability.  The  internal  stability  of  this  loop 
can  be  determined  by  examining  the  four  transfer  functions  shown  in  Figure  B.4. 
These  C2ui  be  written  as 

■  Cl  1  _  f  (/  -  (/  -  AM„)->  A  1  r  u,  1 

.C2  J  “  [(/-M„A)-»Mn  (7-A/„A)-^  J  [  U2  J 

Consider  the  transfer  function  from  U2  to  Ci: 

Cl  =  (/  —  AA/ii)  ^ Au2 

We  can  rewrite  this  as 

Cl  =  AMiiCi  +  Au2. 

Taking  the  2- norm  and  employing  some  properties  of  norms  from  Appendix  A: 

llcilb  =  llAAfiiCi  +  AU2II2  S  ||AAfiiCi|l2  +  IIAU2II2 
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Figure  B.l:  Uncertain  system  with  feedback  control 


Figure  B.2:  Uncertain  system  in  standard  form 


83 


Figure  B.3:  System  expanded  to  show  subsystems  of  Af  (s) 


(B.l) 


lle»lU<llAA/„||oo||e,||2  +  ||A|Ul|u2|h 

||e,||2<l|A|loc||M„|Ul|e,|b  +  llA|UI|t.2lb 

We  can  solve  for  the  2-nonn  of  Ci : 

||e,|U<(l-||A|loc||M„|Uniti2||2 

Noting  that  ||Al|oo  <  li  this  inverse  is  finite  if 

||A/„||oo<l. 

Analysis  on  the  other  transfer  functions  yields  the  same  conclusion.  We  can  now 
state  that  for  the  uncertain  system  of  Figure  B.2  where  the  perturbation  matrix  is 
bounded  such  that 

IIAIloc  <  1, 

internal  stability  is  guaranteed  for  all  allowable  perturbations  provided  that  M{s)  is 
nominally  stable  and 

\\Mu\\oo<l- 
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APPENDIX  C 


Guaranteed  Gain  and  Phase  Margins  Using  the 

Sensitivity  Function 


In  this  appendix  we  derive  one  method  for  providing  guaranteed  gain  and  phase 
margins  for  a  mulit-input  multi-output  system  by  using  a  scalar  input  feedback  un¬ 
certainty  (Dailey,  1990).  Consider  the  plant  in  Figure  C.l.  The  uncertaunty  blocks 
are  arranged  in  a  feedback  loop.  When  the  perturbations  are  combined  into  a 
single  diagonal  block  as  shown  in  Figure  C.2,  the  transfer  function  seen  by  A  is 
(7  —  KG)~^ .  We  recognize  this  as  S,  the  sensitivity  function.  An  equivalent  system 
is  shown  in  Figure  C.3.  In  this  diagram,  A  =  7  —  L  =  KG  '\s  the  loop  gain 
and  C  is  a  diagonal  matrix.  The  gain  seen  by  L  is  equal  to  {1  —  (1  —  =  C. 

Thus,  the  elements  of  C  are  simply  gains  for  each  path  in  the  loop.  From  the  results 
of  the  stability  analysis  in  Appendix  B,  a  guaranteed  upper  bound  on  the  gain  of  A, 
namely  d^(A),  for  the  feedback  loop  to  remain  stable  can  be  determined  from 

d(A)d’(5)  <1  V  w  (C.l) 

Since  the  maximum  singular  value  is  independent  of  the  phase  of  the  matrix  elements, 
e2u:h  diagonal  element  A,  in  Acan  be  expressed  as  a  magnitude  and  phase,  |  A,-  | 
e^,  6  €  [0,27r].  From  Equation  C.l,  d^(A)  <  l/a(S)  is  required  for  stability.  Now, 
solving  for  the  gain  matrix  C  from  above  yields  C  =  (7  —  A)"^  This  matrix  has 
elements  1/(1  —  A,)  on  the  diagonal.  Assuming  only  real  perturbations,  the  gains 
will  vary  from  1/(1-  |  Ai„„  |)  to  1/(1+  |  Ai„.,  1).  Observing  that  for  a  diagonal 
matrix.  A,  o(A)  =1  A,„„  1,  the  stability  requirement  yields 
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Figure  C.l:  Closed  loop  plant  with  input  feedback 


Ai  0 
0  Aj 

5 

Figure  C.2:  Closed  loop  plant  with  diagonal  feedback  block 


Figure  C.3:  Loop  gain 


Thus,  the  gain  margins  aure  given  by  [yi]^ — |  7;^^ — []•  Figure  C.4  shows  a  plot  of 
I  I  versus  gain  mau'gin  in  dB.  It  can  also  be  shown  that  the  guaranteed  lower 

bound  for  the  multi-input  multi-output  phase  margin  is  given  by 

PM  6  \-6,-^d]  where  6  = 

Figure  C.5  shows  a  plot  of  |  |  versus  phase  margin  in  degrees. 
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Figure  C.5:  Phase  margin  in  degrees  versus  | 
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